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(^ ■ Abstract 
(N' 

Vh I We apply a sum rule for the forward light-by- light scattering process within the context of the 

(j)^ quantum field theory. As a consequence of the sum rule a stringent causality criterion is 
presented and the resulting constraints are studied within a particular resummation of graphs. 
Such resummation is demonstrated to be consistent with the sum rule to all orders of pertur- 
. , bation theory. We furthermore show the appearance of particular non-perturbative solutions 
Q_il within such approximation to be a necessary requirement of the sum rule. For a range of values 
A I of the coupling constant, these solutions manifest themselves as a physical bound state and a 
D ' ii'-matrix pole. For another domain however, they appear as tachyon solutions, showing the 
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inconsistency of the approximation in this region. 
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1. Introduction 



Q '. Sum rules provide a powerful tool to study relativistic quantum field theories, and apply 

^^ I also outside the regime where perturbative expansions hold. As sum rules are consequences of 
such general principles as analyticity and unitarity, they allow to establish rigorous relations 
between physical observables, even when the underlying theory is non-perturbative in nature 
rS ■ and cannot be solved exactly. 

In recent works [1, 2], we have derived three sum rules for the low-energy forward light-by- 
light scattering process. These 77-sum rules are non-perturbative in origin and demonstrate 
that the low and high-energy behaviors of the theory are related. We showed e.g. that a sum 
rule for the helicity-difference total cross-section of the photon-photon-fusion process, 77 — !■ X, 
reveals in the hadron sector an intricate correlation between contributions of pseudoscalar and 
tensor mesons. In the charm quark sector, the 77 sum rules reveal an interplay between cc 
bound states, open charmed meson continuum states, as well as exotic cc resonance states [2]. 
Several experiments, primarily at e~^e~ collider facilities, are presently uncovering the rich 
spectroscopy of such systems, see e.g. [3, 4] for some recent reviews. 

Besides its relevance to hadron phenomenology, we can use the 77-sum rules in the same 
way in model field theories, where in the case of renormalizable models they can be applied 
perturbatively. When the conditions of applicability are fulfilled these sum rules were shown 
to hold in leading order calculations [2]. However, the realization of the causality constraints 



Preprint submitted to Physics Letters B April 9, 2013 



at higher orders as well as in the non-perturbative regime of quantum field theories is still an 
open issue. 

Studies of causality constraints on the basis of different sum rules were carried out in the 
past in a number of different contexts. Especially the realization of the well-known Gerasimov- 
Drell-Hearn sum rule [5] within perturbative field theory was analyzed for spin- 1/2 targets at 
the lowest nontrivial order [6] as well as at higher orders in QED [7]. In Refs. [8, 9] consequences 
of the sum rules within asymptotically free theories were considered. In more recent years, they 
have also been discussed within the context of quantum gravity [10, 11]. 

In the present work, we are using light-by-light scattering sum rule as a tool to study causal- 
ity constraints within a model field theory, the 0^ scalar theory. We consider a bubble-chain 
resummation and demonstrate it to be consistent with causality to all orders of perturbation 
theory. Furthermore, it is shown that the sum rule strictly defines the non-perturbative struc- 
ture of the solutions which arise dynamically within this approximation. In a particular regime 
of the coupling constant the spectrum of solutions contains a dynamically generated bound 
state and a i^-matrix pole. For another domain the solution possesses an unphysical pole with 
negative invariant mass being a direct sign of the inconsistency of the approximation. 

The outline of this letter is as follows. In Sect. 2, we compute the light-by-light scattering 
sum rule involving the helicity difference cross section for the 77 — !■ X process, within the 
(j)'^ scalar field theory at one-loop level. In Sect. 3, we provide a calculation beyond the one- 
loop level in the "bubble-chain" approximation. In Sect. 4 we discuss how causality imposes 
constraints on the solutions for different values of the renormalized self-interaction coupling 
constant of the 0^ theory. The summary and outlook are given in Sect. 5. 

2. One loop 

In this work, we will focuss on a sum rule for the forward light-by-light scattering, which 
involves the helicity- difference cross-section for real photons [12, 13, 1] and reads as: 

00 

Acr(s) , , 

ds ^ = 0, (1) 

so 

where Ao"(s) = o"2(s) — cro(s) is the total helicity-difference cross section of the two-photon 
fusion process 77 — )■ X, with the Mandelstam variable s = (g'l + q'2)^, where q\ and g2 are the 
two photon 4-momenta, and Sq is the lowest production threshold of the process. 

We will study the above sum rule in a particular model quantum field theory. We take 
one of the simplest examples: a self-interacting scalar field 0(x) with charge e and mass m as 
described by the following Lagrangian density, 

L = iDyyD,<p - m'<p*<p + ^(0*0)^ - If^'F,,, (2) 

where Aq is the self-interaction coupling constant, while the covariant derivatives and electro- 
magnetic field-strength tensor are given as usual by D^ = d^ + ieA^ and F^j, = d^j^A^, — d^A^. 

We denote the helicity amplitudes for the process 77 — )■ 00* by M4.+ and M+_, where the 
subscripts indicate the photon helicities. Given these amplitudes, the cross section for total 



helicity-0 and 2 are found as: 
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Ms) = ^ / dcose\M+4s,e)\\ (3) 

ozns J_i 



^2(s) = ^ / dcose\M^.{s,e)\\ (4) 

where 6 is the angle of one of the members of the 00* pair w.r.t. the photon in the center-of-mass 
system, and where we introduced their relative velocity /3 as: 



P{s) = ^l-^. (5) 

To leading order in Aq = Ao/(47r)^ and in the fine-structure constant a = e^/47r, the helicity 
amplitudes are found to be: 

MMs,e) = 4vr«{^^^^;|^ + Ao2F(s)}, (6a) 



2/32 sin2 e 



M^.{s,e) = Ana ^ \l _, , . (6b) 



where the first term in M++ and the expression for M+_ correspond with the tree level ampli- 
tudes for 77 — >■ 00* in scalar QED. The second term in M_|__|_ describes the one-loop production 
process corresponding with Fig. 1, with the grey blob denoting the four-particle amplitude 
evaluated to leading order in Aq. Furthermore in Eq. (6a), the (dimensionless) form factor F{s) 
describing the transition of photons to a scalar pair at one-loop order is given explicitely by : 



^2 

F(s) = -1 - ^^Re 
s 



m'„ , 1 + /3(s 



0.-n-/i--i^ ^ 
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It is quite easy now to compute the cross sections for 77 — )■ 00*, the result for the helicity 
difference cross-section is: 

Aa(s) = ^a^'""'\s) -a^\^—[l- p\s)) arctanh/3(s) ReF(s), (8) 

where Ao"'-'^'''"^-' is the tree- level cross section in scalar QED (cf., e.g.. Appendix in [2]). The tree- 
level cross section weighted with 1/s integrates to zero by itself, and it can easily be verified 

that 

00 

ReFis) 
ds ^ arctanh/3(s) = 0. (9) 

Am? 

Hence, we have shown that the sum rule is obeyed at the one-loop level. 



3. Bubble-chain sum 

A class of diagrams naturally arises when one analyzes how the sum rules are realized at 
higher orders in perturbation theory, and in fact beyond perturbation theory. In the following. 






Figure 1: The contribution to the 77-fusion process within the 0^ field theory considered in this work. The 
solid lines denote the charged scalar fields. 





Figure 2: The bubble chain contribution entering the grey blob of Fig. 1. 



we discuss the contribution of the bubble-chain type diagrams to the 77-fusion process as 
shown on Fig. 1, where the shaded blob now denotes a bubble-chain contribution to the four- 
particle vertex, as shown in Fig. 2. The bubble-chain approximation arises in many contexts, 
but most notably as the leading large- A^ result of the 0{N) models. The interest in such an 
approximation is due to the fact that it preserves much of the non-linear structure of the exact 
theory [14, 15]. 

The bubble-chain only contributes to the helicity-0 amplitude M+_|_, with the rath-bubble 
contribution being given by 



r(™+i) 



M|7'^(s) = 47raAo [Ao5(s)]"2F(s) 



(10) 



In the dimensional regularization scheme the one-loop corrections to the four-particle vertex, 
corresponding with a single bubble in Fig. 2, can be expressed through the scalar integral : 



Bis) 



dH 



{2txY [{j> + If -w? + ie] [P -m? + is] ' 
with p the total four-momentum, and s = p^. This integral has the explicit form : 



(11) 



Bis 



-L.ifi^) + 2 - /3{s) In 



i+m 



+ i7iP{s)9{s-4m^) 



(12) 



(47r)2 [ "^^ ^ ' ^^ ^ l-f3{s) 

1/e + 7_E + log(m^/47ryU^) is the dimensional regularization factor and fi is the 



-F'fl) ~ 0.5772 is Euler constant. 



Here L^ = 

corresponding dimensional regularization scale, 7^; 

The resulting amplitudes are independent of the renormalization scale /i which is achieved 
as usual by adding a counter term corresponding to the coupling constant renomalization. We 
define a renormalized coupling constant A by the equation 



A"i(/i2) = V + ^.(/^')-2. 



(13) 



This gap equation allows to renormalize the calculated contribution to each order perturba- 
tively. Using such renormalization, leads to the renormalized (subtracted) one-loop four-point 
function : 



B{s) = (47r)2 [B{s) - B{Am' 



-/3(s)lni±||^ + z7r/3(.)^(s-4m2), 



(14) 



icq 
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Figure 3: Real and imaginary part of the renormalized one-loop correction to the four-point function of Eq. (14), 
as function of s. 



where we have also absorbed a factor (47r)^ in its definition. We notice that our renorniahzation 
procedure is conveniently chosen so as to make a subtraction at threshold: i?(4?Ti^) = 0. We 
show the real and imaginary parts of the renormalized function B{s) versus s in Fig. 3. 

The interference of two chain diagrams with total number of [n — 1) bubble loops gives rise 
to a cross-section correction of the order 0{X^). For the helicity-difference cross-section, which 
in the given case is equal to the helicity-0 cross section, we obtain as correction beyond the 
tree-level : 



Acr(")( 



s 



n-2 



as)Re F{s)B--\s) +|F(.)p^5 



B 



n—2—ii 



(15) 



i=0 



where we used the notation 



as) 



2(l-/3^(^)) 



arctanh/3(s). 



(16) 



One can check explicitly that the expression of Eq. (15) satisfies the helicity-difference sum rule 
exactly in each order of perturbation theory, i.e. 



I(")(A) 



ds 



Aa(")( 



0. 



(17) 



im? 



In Fig. 4, we plot the integrands of the sum rule at different orders of perturbation theory. One 
sees that in all cases the low- and high-energy contributions cancel. 

It is well-known in quantum field theory that when trying to sum the different orders of 
perturbation theory, the resulting series is mostly divergent and has at best a meaning as an 
asymptotic series, a notable exception being the case of theories which have the property of 
asymptotic freedom [16, 17, 18]. As a result, the perturbation theory often can not determine 




Figure 4: The integrand of the sum rule Eq. (17) for the scalar QED with four-particle self-interaction. The 
(red) fine-dashed line corresponds with the one-loop approximation, the (purple and blue) medium- and large- 
dashed lines correspond with two- and three-loop approximations respectively, the (black) solid line corresponds 
with the non-perturbative sum of loop diagrams. The integrands are shown for A~^ ~ 4. 

the solution uniquely. If we are only interested in perturbative phenomena, it is not a problem 
to deal with it in the region of small coupling constants. However if we are interested in non- 
perturbative phenomena then we are faced with the problem to give a meaning to a divergent 
series. 

For example, the photon propagator calculated in the leading logarithmic approximation 
possesses an unphysical pole associated with negative invariant mass and negative probability, 
usually called the "Landau ghost" [19]. In QED, however, due to the smallness of the coupling 
constant the ghost appears at an extremely high energy scale, and the results at the energy 
scales accessible in experiment are not influenced by it. However, in contrast to the QED case, 
in the theory under consideration the coupling constant is not constrained to small values, and 
we are faced with the need to regularize our solution. The appearance of the Landau singularity 
is usually attributed to non-Borel-summability of the considered series, where some individual 
Feynman amplitudes are positive and grow like the factorial of the number of vertices producing 
a singularity on the real positive axis of the Borel transform of the perturbative series [20]. A 
similar situation arises also in the context of our model. We can easily see how non-Borel- 
summability manifests itself in the context of sum rules, by analyzing the contribution from 
the high-energy region to the sum rule integral. In the example under study, the cross section 
Eq. (15) in the ra-bubble approximation behaves at large s as : 

Aa(")(s)~(-A)"^i^^^/^^^^. (18) 



The contribution of the high-energy region to the sum rule integral can then be approximated 
by: 



Act(")( 



-ds 



- (Ing/m^) 

'^) :^- 



-ds 



-A)"(n 



2)!^ 

^ A2 



A2 

In — + 0(l/n) 



VTi^ 



A2 



A2 



(19) 



where A^ ^ vn? . As we see for negative A the high-energy contribution to the sum rule integral 
is positive definite and grows factorially with the order of perturbation theory, which amounts 



to the non-Borel summability of the series 

oo 
oo oo „ 

i(a)=5:i(")(a)=j:/ 



Aor(")( 



ds. 



(20) 



n=0 



n=0 ' 



so 



As a result we can consider ^^q ^'^''"H-^) only in the sense of an asymptotic series. According 
to Poincare [21], a divergent series is an asymptotic expansion of a function /(A) if 



lim 

|A|-i.O 



A^ 



/(A)-5^/,A* 



fc=0 



0, 



for n > 0. 



(21) 



This definition implies that an asymptotic series does not define a function uniquely. The 
expansion coefficients of a function of the type e~^l^ Sfc=o f^^ being identically zero, such a 
function can be added to /(A) without changing Eq. (21)[22]. 

Thus it is natural to expect that in order to obtain the correct behavior of /(A) for negative 
A, one has to modify a formal resummation of the geometric series of renormalized bubble-chain 
corrections to the cross section in Eq. (15), which is given by : 



Aa(s) 



oo 

E 

n=0 



Aa(")(s) = Aa^*'"'^'') - o?—^[s) { e(s)Re 



F{s) 



X- 



B{s) 



F{s) 



X- 



B{s) 



, (22) 



where the tree-level cross section satisfies the sum rule by itself. As was discussed above, due to 
the non-Borel-summability of the series of Eq. (20), we do not have any reasons to expect that 
the sum rule integral will vanish for the resummed theory. In Fig. 5 we show the dependence 
of the sum rule integral for the cross section of Eq. (22) on the value of A. We indeed notice 
from Fig. 5 that the sum rule is only valid for positive values of A (denoted by region I), but is 
violated for negative values of A (regions II and III on Fig. 5), showing that the naive procedure 
of the resummation is not applicable. In order to preserve validity of the sum rules beyond 
the region I, we need to find a way to evaluate the cross section correctly. We will discuss the 
physical situation for the three regions of A in the following. 

0.7 
0.6 
^ 0.5 
0.4 
0.3 
0.2 
0.1 
0.0 i 
-0.1 



b 
<1 





A 




i 


^ 










/ 


<n 




instability / 


-a 


perturbative 
regime 




D 






O 




/ 


X5 




in 


II 


I 


1 


. 





■10 



10 



r' 



Figure 5: The dependence of the sum rule integral for the helicity difference cross section of the 77 
on the inverse coupling A~^. 



X process 



4. Discussion of the results 

Region I : convergent perturbative expansion 

Though the sum of Eq. (20) is formally undetermined, we can still use a naive resummation 
at least in the region of positive A, as one can see from Fig. 5. For A~^ > (region I) the 
series is alternating-sign, since Rei?(s) < 0, and one can expect the series to be resummable. 
Alternatively we can interpret such a resummation at the level of photon-photon-fusion ampli- 
tudes by summing up contributions of bubble-chain diagrams at different orders, which yields 
the amplitude : 

M^4s, 6) = Mf:^\s, 9) + 4vr« ^^^^ , (23) 

with tree- level amplitude given as in Eq. (6a). Squaring the amplitude of Eq. (23) then yields 
the cross section of Eq. (22). In the region I, the amplitude (23) has no poles for all complex 
values of s, and the series is conventionally convergent for all values of s. Thus the formal 
resummation of Eq. (22) satisfies the sum rule. 

Region II : Bound state and K-matrix pole 

Now we proceed to the discussion of the second region, denoted by region II on Fig. 5. In 
the domain of —2 < A"^ < the sum rule is not valid for the naive cross section, given by 
Eq. (22). It is easy to see that the amplitude of Eq. (23) acquires additional singularities below 
the two-particle production threshold. Indeed from Fig. 3 we can see that the imaginary part 
of B{s) vanishes for s < Aim? and monotonically increases for s > Air?. The real part satisfies 
the inequality 

-oo < ReB{s) ^ 0, (24) 

for ^ s < oo, with the upper limit attained for s = 4m^. Thus we see that for A~^ < there 
is always a S'-matrix pole for s < Air?. Above the two-particle threshold, there is a value of s 
where A~^ = Rei?(s), corresponding with a J'i'-matrix pole as will be discussed below. 

We first discuss the pole below the threshold, which obviously corresponds to a bound state 
of the scalar pair. One notices from Fig. 3 that the mass of the bound state varies continuously 
from M^ = Aim? for A~^ = to M^ = for A~^ = —2 as we sweep over the bound state region 
on Fig. 5. In general, the position of the bound state pole is defined by the equation : 

A"^ = B{Ml). (25) 

It is important to note that the bound state singularity is not described by perturbation theory 
being essentially of non-perturbative nature. 




Figure 6: The amplitude of the bound state production in the 77-fusion process. 

The agreement with the sum rule can now be remedied by treating the bound state as a new 
asymptotic state of the theory and thus including the channel of its production in 77 collision, 
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see Fig. 6. The corresponding contribution to the hehcity-0 cross section is: 

aS--^)(.) = 4W^|F(M|)p5(.-M|), (26) 

where the effective couphng of the bound state to 00* is found as the residue of the pole: 

^cff = ~ \ , B\Ml) = ^B{s) . (27) 

'*' \B'{M^)\ ^ "" ds ^ ' s=Ml ^ ' 

It is not difficult to see now that its contribution to the sum rule integral: 

^4--^)(.)=4W|||F(M|)r, (28) 

exactly counter-balances the contribution of the 77 — )■ 00* channel shown in Fig. 5. Thus the 
causality is restored in this region of A. 

We now turn to the position of the singularity above the two-particle threshold. To describe 
the elastic 00 scattering in 0^ theory, we only need to consider S-wave scattering in the bubble- 
chain approximation. The (dimensionless) elastic forward scattering amplitude /(s) is expressed 
through a real phase shift 5{s) as: 

/(s) = e'^(^)sin5(s), (29) 

or through the A'-matrix amplitude K{s) = tan5(s) as: 

/(^) = r,.,\ ■ ■ (30) 

Since the imaginary part of the loop function is 

Im 5(s) = 7r/3(s) {s - Am'^) , (31) 

for s > Am? we can define the elastic amplitude as : 

f{s) = 71(3(3) ^ = (y^ + -arctanh/3(s) - 1] , (32) 

A-i - B{s) \X7rf3{s) vr J 

and hence 

1 2 

K-\s) = -. + -arctanh/3(s). (33) 

A7r/3(s) TT 

In Fig. 7 we show plots of the phase shift for different values of the coupling constant. Note that 
for negative A the phase-shift starts from tt which indicates the presence of one bound state. 
Also for negative A the phase shift crosses 7r/2 at s > Att? satisfying the following equation: 

A-i = Re5(M|.). (34) 

This is the location of the /f-matrix pole, corresponding with a scattering amplitude which 
becomes purely imaginary. Usually this behaviour is attributed to a resonance. Since above the 
threshold, the imaginary part of B{s) is not zero for all complex s and A is real, the amplitude 
(23) does not possess a S'-matrix pole for s > Am?. Hence there is no resonance associated 
with this f^-matrix pole. Note that on the right side of region II on Fig. 5, corresponding 
with A~^ = 0, the f^-matrix pole merges with the bound state and is defined by the position 
Ml- = Air? . When reaching the left side of region II on Fig. 5, corresponding with A~^ = —2, 
the mass is obtained from Rei?(M|;) = 0, which implies M]^ ^ 13.1?n^. 



150 



?? 100 



50 





k 






k 


^^_^:^ 




I 


^^^^^^ --1 




<^ 


■ ____^ A =-2 




"^= 


**===^ r'=o- ~ 


~ 1 






A =0+ 


r 


A"'=16 



10 



20 30 



40 



Figure 7: Phase shift for different values of A. 

Region III : Ground state instability and tachyonic solution 

If A~^ becomes smaller than —2 (corresponding with region III on Fig. 5), the binding energy 
of the bound state exceeds 2?7i, and the pole crosses the point s = 0, moving into the unphysical 
region s < 0, and producing a tachyonic solution (a pole with negative invariant mass). The 
occurrence of a pole for negative values of s signals that the ground state of the theory is 
unstable [14, 15]. We could formally still include this pole in the 77-production channel as an 
asymptotic state and add its contribution to the total cross section in the same way as we did 
for the bound state, which restores the validity of the sum rule. For small negative values of 
A, i.e. A~^ — )■ —00, corresponding with the left asymptotic edge of the region III in Fig. 5, the 
contribution of the tachyon pole is vanishingly small and the sum rule is satisfied approximately. 
This is consistent with the observations in other models, for example, in the leading logarithmic 
approximation in QED, where the position of the Landau pole appears at very large scales due 
to the smallness of the fine structure constant. The position of the tachyon pole in this limit is 
defined by 

-'/\ (35) 



M^ ~ 



which shows explicitly the non-perturbative structure of this contribution. The contribution of 
this pole in this sum rule is asymptotically defined as 



Aai 



1/M^ 



,1/A 



(36) 



and vanishes when A^^ — t- —00. However such a procedure has a number of inconsistencies. 
Being required by the sum rule such extra contribution spoils the general principles of field 
theory. The cross section of Eq. (36) is negative, thus indicating that the appearance of a 
tachyon ghost state contradicts the usual quantum mechanical probabilistic interpretation and 
spoils unitarity. Moreover, the delta function is located at a space-like squared four momentum 
— M^ indicating the occurrence of a tachyonic instability and spoiling the analyticity principle. 
All these facts show that the bubble-chain approximation is essentially inconsistent in region 
III of the coupling constant. For the full theory, such unphysical state will not appear among 
the exact eigenstates. 
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5. Conclusions 

We have studied consequences of causality constraints imposed by a recently established 
sum rule for the forward light-by-light scattering process within the 0^ scalar quantum field 
theory. Within this theory, we verified the sum rule at the one-loop level as well as to all orders 
within the bubble-chain approximation. Furthermore, we have performed a resummation of 
bubble graphs. We have demonstrated that depending on the value of the renormalized self- 
interaction coupling constant of the 0^ theory, three different regimes emerge. In a first regime, 
the perturbative series is convergent and the sum rule as calculated from the continuum states 
in the theory holds exactly. In a second regime for the renormalized coupling, the resummed 
amplitude acquires additional singularities: a dynamically generated bound state below the two- 
particle production threshold and a i^-matrix pole above the two-particle production. It was 
shown that when evaluating the light-by-light sum rule, the bound state contribution exactly 
cancels the continuum contribution, so as to verify the sum rule. Furthermore, we found a 
third regime of the renormalized coupling where a tachyonic solution with negative invariant 
mass appears, signaling that in this regime the vacuum is unstable and that the considered 
bubble-chain resummation is essentially inconsistent. 

The results within the considered model relativistic quantum field theory have demonstrated 
that light-by-light scattering sum rules provide a very powerful tool to quantitatively connect 
dynamically generated bound states with the continuum region of the theory. As such this can 
be a first step, to apply such a tool e.g. to the study of mesons produced in the 77-fusion 
process presently under study at different e"'"e~ collider facilities. 
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